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压 电介质 夹杂 功能 梯度 压 电 带 界面 双 裂 纹 反 平面 问题 " 


周 小 玲 ， 李 Ht 
(宁夏 大 学 数学 计算 机 学 院 ， 银 川 750021) 


摘 要 : 利用 积分 方程 方法 ， 本 文 研究 了 夹 在 两 个 均匀 压 电 半 空间 的 功能 梯度 压 电 带 界面 共 线 双 裂纹 的 
反 平 面 问题 。 在 电 渗透 型 边界 条 件 下 ， 通 过 Fourier 余弦 变换 将 所 考虑 的 问题 化 为 一 对 偶 积分 方 
程 ， 再 用 Copson 方法 将 该 对 偶 积分 方程 转化 为 Fredholm 方程 进行 数值 求解 ， 从 而 给 出 了 裂纹 
尖端 的 应 力 强 度 因子 ， 电 位 移 强度 因子 的 表达 式 。 分 析 了 裂纹 长 度 ， 功 能 梯度 非 均匀 参数 以 及 材 
料 的 几何 尺寸 等 对 应 力 强度 因子 的 影响 。 
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对 于 压 电介质 的 断裂 力学 问题 的 研究 较 充 分 "四 ， 但 是 由 于 压 电 材料 的 本 质 脆 性 ， 导 致 这 
种 材料 的 应 用 受到 了 极 大 的 限制 ， 而 功能 梯度 压 电 材 料 ,"] 较 好 地 克服 了 压 电 材料 的 这 一 次 
端 ， 极 大 地 提高 了 元 件 的 可 靠 性 。 工 程 中 需要 将 不 同 的 功能 梯度 压 电 材 料 拼接 起 来 使 用 ， 所 
以 有 必要 对 此 类 构件 进行 研究 。 然 而 在 拼接 界面 上 极 易 产生 裂纹 并 扩展 ， 此 类 问题 的 研究 较 
多 [B10 。Zhou $19.19 运用 Schmidt 方法 求解 了 两 个 半 无 限 大 压 电 压 磁 复 合 材料 的 界面 含有 两 
个 共 线 界面 裂纹 的 静态 问题 以 及 动态 问题 。 而 运用 Copson 方法 研究 两 种 不 同 的 功能 梯度 压 电 
条 的 界面 含有 共 线 双 裂 纹 拼 接 到 两 个 半 无 限 大 压 电 材 料 上 的 反 平面 问题 还 未 见报 道 


2 问题 描述 


考虑 如 图 1 所 示 的 模型 ， 为 方便 ， 将 功能 梯度 压 电 带 的 上 下 两 部 分 区 域 分 别 用 角 码 1 和 2 表 
示 ， 均 匀 压 电 材 料 上 下 区 域 分 别 用 角 码 3 和 4 表示 。 功 能 梯度 压 电 材料 本 构 方 程 为 


k); OPk e) Owk oth) dpr 

7 = iy) A + P(t, al = cia U) By + els WB (1) 
Ow a ðw 3$ 

DP =a Wa PONa -oy — @) 


其 中 大 = 1,2,3,4, k= 1 2 时 ， 假 设 这 些 材料 系数 沿 y 方 向 按 指数 函数 分 布 ， 即 


ch) (y) = caet, el (y) = eisor, eK (y) = e106, (3) 
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这 里 e440, E110, e150 X y = 0 平面 处 的 压 电 材料 常数 ，Bi 为 梯度 参数 。 而 当天 = 3,4 时 ， 这 些 材 
料 常数 为 


k) 


w 可 +28;h; 
eee elk (y) = el50e Bj i, 


k E 2 
c (y) = case el (y) = eret?” (4) 


其 中 ， 当 k= 3 时 取 正 号 且 7 = 1， 当 k= 4 时 取 负 号 且 7 = 2， 而 hi 和 hs 是 表示 到 有 限 边 界 的 


距离 。 
将 (3)-(4) RRA (1)-(2) 式 ， 得 到 功能 梯度 压 电 材料 的 控制 方程 


V7 we + opp Owe =0, V2wx+ 2 Cu = 0, (5) 
Oy Oy 
XE V? = 07/dx? + 0? /dy? 为 Laplace 算 子 ， 且 有 关系 式 
Wele,y) = delz y) — we (2,y), k= 1,2. 
€110 
均匀 压 电 材料 的 控制 方程 
cd40V2auui + €150V2i = 0， e150 V?wi -elloV2gi = 0, 1 = 3,4. (6) 


Al: 夹 在 两 压 电 半空 间 的 两 不 同 功能 梯度 压 电 带 界 面 双 有 裂纹 的 几何 模型 


3 ” 电 渗 透 型 边界 条 件 的 提出 
(i) 对 称 条 件 
wi(z,0+) = w(x,07), 7 (x,0+) = 72 (z,07), 
$1(x,0+) = $2(2,0-), DP(z,0+) = DP (z,07), lz|>a, lz|<e. (7) 
(ii) 裂纹 面条 件 
Ty2 (7,01) = r (z,07) =—7, fi(z,07) = d(x, 07), 


D(z,0+) = DOGO): c< |z] < a. (8) 
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(iii) 沿 着 界面 的 连续 性 条 件 
TD (a, hp) = T7695, ht), DP(z,h7)=D®(£,h7), w(z, hz) = wz, ht), 
Tya (2,—hd) = {$ (a, ho), D(z, hf) = D(x, —ha*), 
we(x,—hZ) = wa(z,—hz), |z| < œ. (9) 


考虑 到 模型 的 对 称 性 ， 只 需 研 究 z > 0 的 部 分 即 可 。 
Xt FA (8), (9)， 取 关于 z 的 Fourier 余 弦 变 换 ， 再 利用 反 演 变换 得 到 位 移 w(z,y), i = 
1,2,3,4， 与 电势 1 (x,y), i= 1,2,3,4， 的 形式 解 ， 如 


wi(z,y) = a [ Ai(s)e?2¥ + By(s)e"?¥ | cos(sz)ds, (10) 

ġı(z, y) = AG 2) 十 一 sf [Ci(s)e?'y + Dy(s)e??¥ | cos(sx)ds, (11) 
2 fF 入 

w(x, y) = F [A2(s)e™” + B2(s)e*?Y | cos(sx)ds, (12) 
T Jo 


其 中 Ai(s), Bi(s), Ck(s), Dk(s) 是 未 知 函数 ，i = 1,2,3,4; k = 1,2; 
= -fı — V/s? + 3, P2 = -bı + 1/8? + BF, 
Àr = —B2— 9/8? +83, ` = — fo + 1/8? + 63. 


又 由 本 构 方程 (3) 和 (4) BHI r (x,y), i = 1,2,3,4， 和 电位 移 DO (x,y), i =1,2,3,4, 
的 形式 解 ， 如 


Plesy) = Ze 人 [option 十 mmBilsjem) 
+ €150(piC1 (s)e?1¥ + poD; (s)e?2¥) | cos(sx)ds, (13) 


2 œœ l 
D(z, y) = mae / €110 [piCi (s)eP?¥ + poD, (s)e??¥ | cos(sz)ds, (14) 
0 


FEA pp = caio + Seo 


E110 ° 


4 “对 偶 积分 方程 的 建立 
取 裂 纹 上 下 表面 的 位 移 差 和 电势 差 为 
f(x) = wi(x,0*) — we(2,07), F(z) = dr(x,0*) — $2(z,07). (15) 
将 位 移 ， 电 势 的 形式 解 (10), (11) 等 代入 方程 (15)， 并 经 Fourier 余 弦 变 换 ， 可 以 得 到 
f(s) = Ai(s) + Bi(s) — Aa(s) — Ba(s), 


(s) = 252 F(s) + C1(s) + Di(s) — C2(s) — Da(s) = 0, (16) 
E110 


Sl 
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应 用 边界 条 件 (7)-(9)， 经 过 Fourier 余弦 变换 ， 从 位 移 ， 电 势 的 形式 解 (10), (11) 等 和 应 力 ， 
电位 移 的 形式 解 (13), (14) 等 可 以 得 到 到 一 组 关于 A(s), Bk(s), Ck(s), Dk(s) 的 方程 ， 并 与 方 


程 (16) 联 立 可 以 得 到 
Ax(s) = qi(s)q(s)f(s), 42(s) = 92(s)aa(s)a(s) f(s); 
As(s) = [qi(s)e® + + elt)" a(s) f(s), 
Aa(s) = [qa(s)e(°"2 + ef 2)" | qa(s)a(s) f(s), 
Bi(s) = 4(s)f(s), B2(s) = a3(s)a(s) f(s), 


e x — 
B3(s) = e la (s)e(Pits)ha 十 elP2+ )hi ]e(s)7(s)， 


Ba(s) = ae [a2(s)e(s>h2 + el->a)ma | g5(s)q(s)F(s), 


Ci(s) = -2 a(o) FC), C2(s) = go(s)q3(s)a(s) f(s), 
Dy(s) = ~ ~ Sw s)F(s); i 
其 中 gi(s), i = 12,3， 均 是 已 知 函 数 ， 分 别 为 


n(s) = pı (s) +s 


q(s) = 1/[q1(s) + 1 — a2(s)qa(s) — qa(s)]. 
由 裂纹 面 处 受 反 平面 剪 切 力 ， 上 下 表面 应 力 相 等 ， 得 对 偶 积 分 方程 


TO 


Jo” sa(s) f(s) cos(sx)ds = 一 TA <a, 
Jo” f(s) cos(sx)ds = 0, r>a, O<T<C, 
其 中 af(s) 是 已 知 函 数 ， 且 
a(s) = -E + p(s) = PENCE) + e200)las), Jim a(s) = -3. 


通过 简单 分 式 变换 x = z* + ES, b = 号 < 为 方便 略 去 * 后 方程 (27) BAH 


0 


| Jo. sa(s)f(s) cos(sx)ds = oor Iz| < b, 


Jo” F(s) cos(sz)ds = 0, |z| >b. 


对 于 对 偶 积 分 方程 (29)， 我 们 采用 Copsonl A 方法 求解 ， 分 情况 进行 讨论 。 
情况 1 当 z > 0 时 ， 引 入 函数 pf6)， 使 


b 
I(s) = f ép(€) Jo(sé)de, 


agi, gga ileo, a- maara 


(17) 
(18) 
(19) 
(20) 


(21) 
(22) 
(23) 


(24) 


(25) 


(26) 


(27) 


(28) 


(29) 


(30) 
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其 中 jh( ) 为 第 一 类 0 阶 Bessel 函数 。 将 式 (30) 代入 (29)， 则 第 二 个 方程 已 经 自动 满足 ， 要 使 
第 一 个 方程 也 成 立 ， 则 8(&) 必须 由 下 面 的 第 二 类 Fredholm 积分 方程 控制 ， 即 


v(é) + [ CR = -j OS BSS, (31) 
其 中 的 核 函 数 KE, n) A 
Klem) =n | sials) ~ (on) Jo(sé)as. (32) 
为 了 便于 数值 求解 ， 再 次 引入 下 面 的 无 量 纲 变量 
o= 2, €=b2, n=bl, z= Xb, (33) 
oo- on) =- 22 22, (34) 


将 式 (33), (34), RADZE (31), (32)， 则 各 项 分 别 标准 化 为 


1 
FXS) = Bud Í VGW(2) Jo(SQ)aQ, (35) 
W9) + f Fo Dear = VA, (36) 
其 中 
LQ, T) = VAT J 7 s[a(?) — 1] Jo(ST)Jo(S)28, (37) 
0 
a(2) _ al +6(2) = ei pate) (38) 


情况 2 对 于 求解 zx < 0 的 情况 与 情况 1 中 的 做 法 一 致 ， 故 此 处 省 略 。 


5 “裂纹 尖端 电 弹 性 奇异 场 


对 于 断裂 力学 来 说 ， 我 们 感 兴趣 的 是 裂纹 尖端 附近 的 应 力 ， 应 力 强度 因子 ， 电 位 移 以 及 电 
位 移 强度 因子 。 因 此 沿 着 裂纹 面 的 应 力 Ty 和 电位 移 Dy BX > 1 时 可 表示 为 


Tyz(z,0) = cao sa(s) f(s) cos(sx)ds 
1 o0 
= mg f Jh(S) coe(SX)dS +--+, (39) 


Dy, (x, 0) = “iso 1 sa(s)f(s) cos(sx)ds 


= £150 mg() f 1(S) cos(SX)dS +--+ , (40) 
C440 0 
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当 关 < 一 1 时 有 类 似 表 达 式 ， 这 里 (..) 是 第 一 类 一 阶 的 Bessel 函数 ， 亚 (1), U(-1) PAA W(Q) 
在 裂纹 左右 尖端 Q = 土 1 处 的 值 。 由 于 方程 (39), (40) 等 中 的 积分 在 z = 二 处 是 有 界 的 ， 
此 ， 应 力 和 电位 移 的 奇异 性 只 由 包含 亚 (1), Y-I 的 项 控制 。 记 左右 两 端 应 力 强 度 因子 分 别 
为 KT 和 KE， 电 位 移 强 度 因 子 KP 和 KR， 按 通常 应 力 强度 因子 的 定义 可 得 


K= lim V2n(ce—2)tyz(2,0), KR= lim, 2n(z ~ a)Tyz(£,0), (41) 
KP = lim V2n(c~2)D,(z,0), KR = lim V2r(z ~ a)D,(z,0), (42) 
则 它们 可 分 别 表示 为 
1 
KE = -VT im W(-1), KR= 3Vibroy(), (43) 


KP = -£ Vib img h(—1) = 22 Kt, KR = Vrbroy(-1) = FE KR, (44) 
2c440 C440 2c440 C440 


其 中 亚 (1) 与 于 (~1) 分 别 为 更 (O) 在 裂纹 尖端 Q = 1 的 值 。 标 准 化 应 力 强度 因子 


KF 1 
Ke = = = -3iV(- 1), KR = ek = 540). (45) 
标准 化 电位 移 强度 因子 
Kp’ = Ke. £180 g(a), pt = KE = £90 yay, (9 


VrbDy 2c440 D VrbDo 2c440 Do 


由 式 (43), (44)， 可 以 看 出 ， 当 9 一 土 1 时 ， 即 无 限 接近 裂纹 尖端 时 ， 无 限 长 功能 梯度 压 电 材料 
界面 共 线 双 裂纹 尖端 附近 的 应 力 ， 电 位 移 都 具有 才 阶 的 奇异 性 。 


6 ”数值 算 例 与 结论 


BE y = 0 平面 处 的 材料 为 压 电 陶 瓷 PZT-5H， 其 材料 常数 分 别 为 cd4 = 3.53 x 1010 Nm-2， 
el5 = 17.0Cm™, e11 = 150.4 x 10-10C/Vm， 假 定 应 力 载荷 加 = 1.2 x 105 N/m? 及 电位 移 载 
fii Do = 0.01 C/m”。 由 于 电位 移 强 度 因 子 与 应 力 强度 因子 存在 线性 关系 ， 所 以 我 们 只 给 出 应 力 
强度 因子 随 材 料 的 几何 尺寸 和 材料 的 非 均匀 参数 的 变化 情况 ， 有 关 结 果 可 用 图 2 至 图 5 表示 。 


7 ”结论 


分 析 图 2 至 图 5， 我 们 可 以 得 到 如 下 结论 : 

1) 应 力 强度 因子 与 裂纹 间距 ， 材 料 非 均匀 参数 ， 材 料 几 何 尺 十 有 关 ， 但 与 材料 本 身 的 常 
数 无 关 。 

2) 应 力 强度 因子 Ks 随 着 裂纹 半 长 度 b 的 增 大 而 减 小 。 

3) 应 力 强 度 因 子 Ks 随 着 上 层 功 能 梯度 压 电 带 宽度 hl 的 变 大 而 减 小 ， 逐 步 趋 于 一 个 常 
值 ， 市 随 着 下 层 功 能 梯度 压 电 带 宽度 ho 的 增 大 而 增 大 ， 逐 步 趋 于 一 个 常 值 ， 因 此 适度 控制 材 
料 宽度 可 以 减 小 应 力 强 度 因 子 。 

4) 应 力 强 度 因 子 随 着 材料 的 非 均匀 参数 81 的 变 大 而 逐渐 变 大 ， 而 随 着 Bo 的 变 大 而 减 小 ， 
因此 调节 材料 的 非 均 匀 参 数 对 应 力 强度 因子 的 影响 也 是 显著 的 。 
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图 2: 标准 化 应 力 强度 因子 Ks 随 Bi (b RE) 和 随 Bi (hy TE) 的 变化 情况 


5y r x 一 3.5 < — — 一 一 
4.5} 
7 . a 日 口 a b 
了 g 3 = a 
人 
中 baer á o 
3.5F o -Q È ; esto: o o © a RE T, ? 
o o 2 io gee 
ab OO @ pare eres o- Ka B=10b | ge 
mp. 
ca EA —a— Ke T 
Ks 2 5b a a is} s] a 2 O- K 2 o. Ka BB 4 
2} ~O— Ke 
ace 1.5 站 Ka B=0.5B 4 
oeeo o o eo —$— ke 
G_6-__9-—6 - e 
5—6 ~ 全 3 所 一 
i ais 


sas Ka Be=10B, 


一 中 一 Ke 
一 一 - 


图 4: 标准 化 应 力 强度 因子 Ka Ri hy (b 不 同 ) Ri hi (B2 不 同 ) 的 变化 情况 
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Ka 


15 ren 
1.4 上 
1.37 日 Kahi=0.2,h2=0.2 
| —a— ke 
12 o- Kaht=1,he=1 
11 § | —o— Ke 
g- Kahi=5,he=6 


0.5 0.6 0.7 0.8 0.9 1 


图 5: 标准 化 应 力 强度 因子 Ka 随 6 (hi, ho 不 同 ) OM b (Go 不同 ) 的 变化 情况 
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The Anti-plane Collinear Cracks Problem of Functionally Graded 
Piezoelectric Strips Bonded to Homogeneous Piezoelectric Media 


ZHOU Xiao-ling, LI Xing 


(School of Mathematics and Computer Sciences, Ningxia University, Yinchuan 750021) 


Abstract: By using the integral equation method, we investigate the anti-plane problem of two 
collinear cracks between two different functionally graded piezoelectric strips, sandwiched by two ho- 
mogeneous piezoelectric media under shear load. By using the Fourier cosine transformation, the 
considered problem could be reduced to a pair of dual integral equations which are numerically solved 
by the Copson method under the permeable boundary conditions. The expressions for the stress inten- 
sity factor and the electric displacement intensity factor are obtained. Numerical examples show that 
the effects of the length of crack and material gradient parameter, etc. 

Keywords: functionally graded piezoelectric strip; permeable crack; Fourier cosine transform; Copson 


method; dual integral equation; stress intensity factor 
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